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ABSTRACT 

In  this  paper,  we  consider  the  random  effect  panel  data  model  which  has  fixed  and  random  effects  as  well  as 
the  experimental  error  term.  Bayesian  approach  employed  to  make  inferences  on  the  model  coefficients.  To  illustrated  the 
effectiveness  of  the  methodology.  We  have  chosen  a  data  set  from  gross  fixed  capital  formation  and  gross  domestic 
product  by  economic  activities  for  public  sector  at  current  prices  for  the  years  (2005-2015)  (Million  I.D.).The  data  are 
analysed  according  to  our  methodology  by  using  gretl,  R  and  matlab  softwares. 
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INTRODUCTION 

Linear  models  play  a  central  part  in  modern  statistical  methods.  On  the  one  hand,  these  models  are  able  to 
approximate  a  large  amount  of  metric  data  structures  in  their  entire  range  of  definition  or  at  least  piecewise.  The  theory  of 
generalized  models  enables  us,  through  appropriate  link  function,  to  apprehend  error  structures  that  deviate  from  the 
normal  distribution,  hence  ensuring  that  a  linear  model  is  maintained  in  principle.  Linear  statistical  methods  are  widely 
used  as  part  of  this  learning  process.  In  the  biological,  physical,  and  social  sciences,  as  well  as  in  business  and  engineering, 
linear  models  are  used  in  both  the  planning  stages  of  research  and  analysis  of  the  resulting  data.  As  well  as  and  to  the  best 
for  our  knowledge  the  Linear  models  and  Bayesian  models  were  studied  by  many  researchers  for  example 
see, [2], [3], [4], [5], [7], [8], [11], [12], [13], [14], [16]. 

In  this  paper,  we  consider  the  random  effect  panel  data  model.  Our  paper  is  related  to  the  previous  works  [9],  [10], 
which  provides  of  theoretical  results  for  panel  data  model  as  well  as  Bayes  panel  data  model  based  on  Markov  Chain 
Monte  Carlo  (MCMC).  Consider  the  model 

Yit  =  p  +  ZjLx  PjXjit  +  eit, i  =  1 . N,t  =  1 . T,  (1) 

Where,  Yit  the  value  of  response  variable  for  ith  unit  at  time  t,  Xjitthe  explanatory  variables,  p, p.,j  =  1,  ...,K  are  fixed 
parameters  and  eit  is  an  error  term  with  sitn5N(0,  a^). 

Now,  if  the  parameter  p  is  specified  as: 

li  =  P  o  +  uu  (2) 

Where,  Ui~N(0,  o^)t  then,  the  model  (1)  is 
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Yit  =  P0+£f=i  Pjxjit  +  ui  +  £it-  (3) 

The  model  (3)  is  rewrite  as  follows 

Yit  =  P0  +  i  Pjxjit  +  Mitt  (4) 

Where, coit  =  ut  +  sit,coit~N(0,  of  ),<jf  =  +  of,  thus  by  using  matrix  notation  the  model  (4)  is 

Y  =  F0  +  co  (5) 

Where,  F  =  [e,X\,e  =  [  1,1,  has  length  NTY  =  [Yllf ... ,  Y1T,  Y21, ... , ...  Y2T, ... ,  YN1, ... ,  YNT ]T  has  length  NT, 

X=  [Xx,X2,  ...,XN  ]T  is  a  NT  x  K  design  matrix  of  fixed  effects,  6  =  [PQ,  ,  PK]T  has  length  K  +  1,  and  co  = 

[c011;  ... ,  C01T,  (021,  ■■■  ®2T>  ■■■  >  ®N1'  ■■■ '  ®NT  ] 

Has  length  NT.  From  model  (5),  we  haveY~  N  (  F  6,  V  ),  where 
^  =  E  {coco7)  =  IN  ®(o£2  It  +  ol  eeT)  =  ol(IN®It)  +  o2{IN  ®eeT ), 

Replace  Ixby  (ET  +  JT)  and  eeT  by  T  JT,  where  JT  =^eeT  and  ET  =  IT  —  JT,  then 
y=  o2£[IN®{ET+JT)]  +  o2(In®T  Jt) 

=  o2{In®Et  )  +  g£(Jn®Jt)  +  To2{In®Jt ), 

By  collecting  terms  with  the  same  matrices,  we  get 

¥  =  <t£(In®Et  )  +  (a2  +  Tct2)(In®Jt  )  =  o2  Q  +  ol  P,  where, of  =  (o£2  +  Tof)  and 

XF~1  =  I'TI  =  product  of  its  characteristic  roots,  [1]  ->  l^l  =  {ol)N(<T~1>}{ol)N . 

°£  G1 

The  likelihood  function  is  the  joint  density  of  the  Y's  that  is 

—NT  -1 

L(Y-,e,W)  =  (2n)—\V\~T exp{-±(Y  -  Fdyv-'iY  -  Fd)} 

-NT  -N(T- 1)  —N  r  n  Pi  , 

=  (2? r)  2  (cr|)  2  (cjZ)  2  exp  {  —  (Y-  FO  )T  +  -j  (7-  FO  )}. 

Then,  the  likelihood  estimators  of  parameters  0,g2,  a2  are  [9] 

e  =  (Fr¥_1  F  )_1  (  F^-1  7  ),  <t|  =  (7  -  F  0)T  Q  (Y  -  F§)  and  a?  =  ±  (Y  -  F  df  P  (  (7  -  F§) 

The  panel  data  model  has  been  investigated  by  many  researchers  for  example  see  [1],[6  ],[9  ],  [10  ],  [15  ]. 
To  illustrated  the  effectiveness  of  the  methodology.  We  have  chosen  a  data  set  from  gross  fixed  capital  formation  and  gross 
domestic  product  by  economic  activities  for  public  sector  at  current  prices  for  the  years  (2005-2015)  (Million  I.D.).The 
data  are  analysed  according  to  our  methodology  by  using  gretl,  R  and  matlab  softwares. 

The  Prior  and  Posterior  Distributions 

To  specify  a  complete  Bayesian  model,  we  need  a  prior  distribution  on (#,  of,  cr2).  We  will  use  the  uniform 
distribution  U(0,1)  of  the  vector  parameters#,  as  well  as  we  will  assume  that  the  prior  distribution  on  ol  and  of  are  invers 
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gamma  with  of  parameters  a£ ,  pg,  o^and  respectively,  [10] 

exp(-|)  and 


n0  (ft)=:?£:(ftr(ae+1) 


r(aE) 


n0(a?)  =  ^  (ft)  exp  (  —  f),  where  ae  ,p  ,  a.! and  (3  are  hyperparameters  that  determine  the  priors  and  must  be 

FyCt  i)  V  o^J  8  1 

chosen  by  the  statistician. 

From  the  model  (5)  we  have  Y\6,  cr|,  F6,  V  ).  Then,  the  likelihood  function  is 

L(F|0,ft  ,ft)  =  |<P|ftexp{-;(F  -  F0)r'F“1(Y  -  F0)}. 

t=  1 

In  the  exponent,  we  add  and  subtract  F§  to  obtain,  [10] 

[(F  -  Fey  IF_1(F  -  F0)]  =  [(F  -F9  +  F9-  F0)r'F~1(F  -  F0  +  F9  -  Fd )] 

=  [(f  -  Fe)  -  f(6>  -  S^V^F  -  f 0)  -  f(0  -  0 )] 

=  [(Y  -  F0)T'F_1(Y  -  F e)  -  (Y  -  F0)r'F“1F(0  -  9)  - 
{9  -  0)rFT'F“1(Y  -  F0)  +  (0  -  0)rFT'P_1F(0  -  0)], 

Since,  (Ft'F“1  F)0  =  FT'F_1  Y,  then 

[(Y  -  Fe)T'P“1(Y  -  F0)]  =  (Y  -  F0)TvF_1(Y  -  F0)  +  (0  -  §)TFTxP~1F(d  -  0). 

The  joint  posterior  density  of  the  coefficients  0  and  the  variances  ft  andft  givenby  the  expression 
7Ii(0<°e>°'l  IO  K  1(F|0,  ft ,  ft)7To(0,  ft ,  ft) 


—NT  -N(T- 1) 

oc  (2n)  2  (<j£2)  2 


-N  1  T  (  0  P  \  1  T 

(ft)”  exp{—  -  (F  -  F 0)  +  —  (F  -  F0)}  exp  {-  -  (0  -  0)  Fr<F_1F 

z  \oE  ox )  z 


ft 


If) 


-ft 


exv\—\x  (ft)  (“1+1) exP  { -ff ) 


* (ft) 


-(a£+^+1) 


2  +1)  exp  —  ’ 2 


\  (Y  -  Fey Q{Y  -  F9)  +  PE 


x  (0-2) 


-(a1+f+1) 


(-  (F  -  F9)tP(Y  -  FG)  +  ft)  1  ^  T 

exp-j- - 2 - [exp  {--(0  —  0)  FtT  XF(0  —  0)}. 

I  CTi  I  Z 


From  this  expression,  we  can  deduce  the  following  conditional  and  marginal  posterior  distributions 
jri(0|ft,o2,Y)a  exp  {-  ±(0  -  0)TFtT-1F(0  -  0)}, 

And 


(6) 


(7) 
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ni(c£\9,(j(,Y)  oc  (ft 


,  2,a  2  ^  y  2^(ai+-+l)  f  5(i'-fe)TP(y-Fe)  +  /?1'| 

ttiOftlftft.F)  oc  (ft)  v  2  ftxpj-2 - ^ - j 


Therefore,  it  follows  that 

( e\alalY)~N{6 ,  (Fr<F_1F)_1), 

(<t£2  |  0,  a2,  F)~/G  (a£  +  Be  +\(Y-  F9f  Q  (Y  -  F9  )), 

(ft|ftft,F)~/G  (a,  +  f,ft  +  ^  (F  -  F9fP  (F  -  Fdj) 

Bayes  Factor 

We  would  like  to  choose  between  a  fully  Bayesian  panel  data  model  with  ( K  +  l)of  parameters  against  a 
Bayesian  panel  data  model  with  ( q  +  l)ofparameters,  where  q  <  K,  by  using  Bayes  factor  for  two  hypotheses 


(8) 

(9) 

(10) 
(11) 

(12) 


co 


Ho-  Yit  =  Bo  +  Tj=1PjXjit  +  coit,or  H0 :  F°9°  + 
against  ' 

ft: Yit  =  Bo  +  Ey=i PjXju  +  colt, or  ft:  F0  +  co  J 


(13) 


Where,  6°  is  (<7  +  1)  vectors  of  parameters,  F°  is  an  NT  X  (q  +  1)  design  matrix  and  q  <  K.  We  compute  the  Bayes 
factor,  £0i  of  H0  relative  to  Ht  for  testing  problem  (13)  as  follows 


— 


m(Y\H0) 


01  mWi)’ 

Where  m  (Y\Hi)  is  the  marginal  density  of  Funder  model  Hifi  =  0,1.  From  [10]  we  have: 

m(Y\Ho)  =  JJ  (J  f  (Y\e°,oi,o^n1(e0\oi,ol)n0(oi,ol)de°)doidol). 

-NT  ft*1  ft“£  (N(T  - 1)  \  /I  \ 

=  (2")^ffe)7=fer(-4-!  +  “*  +  2)(z  (i" - F°eVW - +  a) 

if)  +  «  +  2)  f*  (y  -  f°9°)tp(y  -  f°9°)  +  ft)  (l,l+,+11 

m(F|ft)  =  //(//  (Fft,ft,ft)7rftftft,ft)7ro(ft,ft)d0)dftdft. 

-NT  ft*1  Be*  /H(T  - 1)  \  rl 

=  (27r)^F(feFfcr  (—2—  +  a*+2)  [2  -  F6?)W  -  Fe)  +  A 

X  r  g  +  a,  +  2)  (i  (F  -  F9)tP(Y  -  F9 )  + 


(14) 


-(a, 


^^0 


_(a  +^0^12+!) 

Q  (F  -  F°9°)tQ(Y  -  F°9° )  +  ft)  V  £  2  '  Q  (F  -  F°9°)tP(Y  -  F°9° )  +  ft) 


B, 


-(a1+f+l) 


01 


Q  (F  -  F9)TQ(Y  -  F0)  +  Be)  2  Q  (F  -  F0ftP(F  -  Ffl)  +  ft) 


-(«1+f+r) 
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The  Data  Results 

To  illustrated  the  effectiveness  of  the  methodology.  We  have  chosen  a  data  set  from  gross  fixed  capital  formation 
and  gross  domestic  product  by  economic  activities  for  public  sector  at  current  prices  for  the  years  (2005-2015)  (Million 
I.D.).The  data  are  analysed  according  to  our  methodology  by  using  gretl,  R  and  matlab  softwares.  We  used  gross  fixed 
capital  formation  at  current  prices  for  the  years  (2005-2015)  (Million  I.D.)  for  all  economic  activities  as  a  dependent 
variables,  gross  fixed  capital  formation  at  current  prices  for  the  years  (2005-2015)  (Million  I.D.)  for  previous  year  and 
gross  domestic  product  by  economic  activities  for  public  sector  at  current  prices  for  the  years  (2005-2015)  (Million  I.D.)  as 
explanatory  variables,  as  well  as,  we  consider  the  economic  activities  as  sections.  Then  we  have  (9)  sections,  where 

Section  (1):  agriculture,  forestry,  hunting  and  fishing, 

Section  (2):  mining  and  quarrying, 

Section  (3):  manufacturing  industry, 

Section  (4):  electricity  and  water, 

Section  (5):  building  and  construction, 

Section  (6):  transport  communications  and  storage, 

Section  (7):  wholesale,  retail  trade,  hotels  and  others, 

Section  (8):  banks  and  insurance, 

Section  (9):  social  and  personal  services. 

This  application  has  been  divided  into  two  parts.  The  first  part  includes  the  estimation  of  sub-models  for  sections. 
By  using  gretl  software  we  obtain  the  maximum  likelihood  estimators  for  the  sub-models  of  the  sections  both  individually 
and  the  total  model  for  all  sections  together.  Table  (1)  below  shows  the  results  for  the  model  estimators.  From  this  table  we 
can  see  that  the  models  for  building  and  construction  and  wholesale,  retail  trade,  hotels  and  others  where  significance  at  a 
0.05  level  of  significance.  This  means  there  is  a  significant  effect  for  gross  fixed  capital  formation  at  current  prices  for  the 
years  (2005-2015)  (Million  I.D.)  for  previous  year  and  gross  domestic  product  by  economic  activities  for  public  sector  at 
current  prices  for  the  years  (2005-2015)  (Million  I.D.)  on  the  gross  fixed  capital  formation  at  current  prices  for  the  years 
(2005-2015)  (Million  I.D.)  for  all  economic  activities.  Furthermore,  for  the  sections(l),  (2), (3),  (4), (6), (8), (9)  which  was 
non-significant,  in  fact  this  don't  agree  with  the  economic  theory,  then  we  can  treat  this  problem  by  using  confounding 
approach  for  this  sections  data  (i.e.  the  data  of  the  economic  activities:  agriculture,  forestry,  hunting  and  fishing,  mining 
and  quarrying,  manufacturing  industry,  electricity  and  water,  transport,  communication  and  storage,  banks  and  insurance, 
social  and  personal  services)  with  time  series  to  obtain  (29)  observation.  Table  (2)  shows  the  results  the  total  model 
estimator  for  all  the  economic  activities,  clearly  from  this  table  the  model  is  significance  at  a  0.05  level  of  significance 
according  to  F-  value  which  (20.27009)  with  p-  value  (F)  which  (7.00e-19),  as  well  as  the  value  of  determinant  coefficient 
( R 2  =  0.70),  this  means  the  total  model  is  agree  with  the  economic  theory.  Also  table  (3)  presents  summary  values 
estimation  of  the  confounded  model.  From  this  table  we  can  see  that  the  total  model  after  confounding  also  significance  at 
a  0.05  level  of  significance.  The  second  part  we  applied  our  methodology  (Bayesian  method  )  to  the  data  set  from  gross 
fixed  capital  formation  and  gross  domestic  product  by  economic  activities  for  public  sector  at  current  prices  for  the  years 
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(2005-2015)  (Million  I.D.).  Figure(l)  represents  the  posterior  density  of  the  coefficients  for  the  model  of  section(l) 
(building  and  construction)  and  figure  (2)  shows  the  number  of  iterations  of  the  Gibbs  sampler  used  in  this  study,  which 
was  (7000  )  iterations,  while  figure  (3)  shows  density  estimates  based  on  (7000  )  iterations  of  a2, a  2  and  for  this  model. 
We  have  displayed  posterior  density  of  coefficients  for  the  model  of  section  (1)  (wholesale,  retail  trade,  hotel  and  others)  in 
figure(4)  and  figure  (5)  shows  of  iterations  of  the  Gibbs  sampler  of  this  model  which  was  (7000  ).  Figure  (6)  represents  the 
posterior  density  based  on  (7000)  iterations  of  a2,  a  2  and  a £  for  the  model  of  section  (2).  Furthermore,  the  figures  (7)  and 
(10)  represents  the  posterior  density  of  the  coefficients  for  the  total  model  and  the  confounded  model  respectively,  also  the 
figures  (8)  and  (11)  shows  the  number  for  the  iterations  of  the  Gibbs  sampler  which  was  (5000)  and  (8000)  for  this  models 
respectively  and  figures  (9)  and  (12)  shows  the  posterior  density  based  on  (5000  )and(8000)  iterations  of  a2,  a2  and  for 
the  total  model  and  the  confounded  model  respectively.  Table  (4)  presents  the  values  of  the  parameters  for  the  models 
estimates  based  on  the  Bayesian  method.  From  this  table,  we  can  see  that  the  values  of  parameters  obtained  by  Bayesian 
method  are  encouraging. 


Table  1:  The  Model  Estimators  for  Sections  (Economic  Activities) 


Section 

Parameter 

estimator 

The  value  of 
estimators 

Standard 

error 

t-  test 

P- 

value 

F  -  test 

p  —value 

(F) 

R2 

1 

0o 

3.86660e+07 

2.55386e+08 

0.1514 

0.8834 

0.213435 

0.812260 

0.050656 

01 

0.0675264 

0.335536 

0.2012 

0.8455 

P2 

17.0406 

28.9661 

0.5883 

0.5726 

2 

0o 

-2.14066e+08 

5.75331e+08 

0.3721 

0.7195 

0.785175 

0.488257 

0.164085 

Pi 

0.0835356 

0.163748 

0.5101 

0.6237 

02 

12.2181 

9.76379 

1.251 

0.2462 

3 

Po 

1.10795e+010 

6.58628e+09 

1.682 

0.1310 

0.880369 

0.451263 

0.180390 

Pi 

-2.45741 

3.28918 

-0.7471 

0.4764 

P2 

-  1324.87 

1456.62 

-0.9095 

0.3896 

4 

Po 

1.28578e+09 

1.52223e+09 

0.8447 

0.4228 

4.093918 

0.059649 

0.505802 

Pi 

0.606713 

0.237005 

2.560 

0.0337 

P2 

170.653 

137.960 

1.237 

0.2512 

5 

Po 

-2.00648e+08 

1.60905e+08 

-1.247 

0.2477 

7.638054 

0.013955 

0.656300 

Pi 

-0.494557 

0.301930 

-1.638 

0.1401 

Pi 

64.7116 

17.6098 

3.675 

0.0063 

6 

Po 

6.51350e+08 

1.64660e+09 

0.3956 

0.7028 

1.215018 

0.346112 

0.232984 

Pi 

0.386163 

0.504278 

0.7658 

0.4658 

02 

37.8384 

181.875 

0.2080 

0.8404 

7 

00 

2.02903e+08 

7.13551e+07 

2.844 

0.0217 

12.04065 

0.003867 

0.750634 

0i 

1.18771 

0.261094 

4.549 

0.0019 

02 

-18.9722 

7.96265 

2.383 

0.0444 

8 

0o 

9.36255e+07 

2.92096e+07 

3.205 

0.0125 

2.066738 

0.188981 

0.340667 

0i 

-0.790978 

0.404201 

-1.957 

0.0861 

02 

12.5969 

9.70421 

1.298 

0.2304 

9 

00 

6.70564e+09 

5.63077e+09 

1.191 

0.2678 

1.239066 

0.339801 

0.236505 

0i 

0.420409 

0.301515 

1.394 

0.2007 

02 

52.6157 

197.660 

0.2662 

0.7968 
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Table  2:  The  Total  Model  Estimator 


Parameter 

Estimator 

The  Value  of 
Estimators 

Standard 

Deviation 

t- 

Test 

p  —Value 

F  —Test 

p  — Value(F) 

R2 

A> 

8.26169e+08 

5.45353e+08 

1.515 

0.1334 

k 

0.390519 

0.0933341 

4.184 

6.76e-05 

20.27009 

7.00e-19 

0.697283 

k 

48.4406 

26.5657 

1.823 

0.0716 

Table  3:  The  Total  Model  Estimator  After  Confounded 


Parameter 

Estimator 

The  Value  of 
Estimators 

Standard 

Deviation 

t  —Test 

p  —Value 

F  —Test 

p  — Value(F) 

R2 

A, 

-1.21363e+08 

9.2781  le+07 

-1.308 

0.2023 

k 

1.04514 

0.00878075 

119.0 

4.06e-037 

9505.042 

5.76e-38 

0.998634 

k 

6.17952 

0.540683 

11.43 

1.23e-011 

Table  4:  Estimation  Values  by  Bayesian  Method  for  the  Models 


The  Model 

The  Parameter 

The  Value  of  Parameter 

Model  of  section  1 

A) 

-2.00626e+08 

k 

-0.494 

k 

64.71 

Model  of  section2 

A> 

2.02091e+08 

k 

1.19 

k 

-18.97 

Total  model 

Ai 

8262e+08 

k 

0.3905 

k 

48.44 

Confounded  model 

A> 

-1214e+08 

k 

1.045 

k 

6.18 

The  model  checking  approach  based  on  Bayes  factors  has  been  tested  on  estimated  models.  These  Bayes  factors 
are  given  in  table  (5).  From  this  table,  it  can  be  seen  that  the  Bayes  factors  favors  H1  with  strong  evidence  with  all  models 
for  the  data  of  gross  fixed  capital  formation  and  gross  domestic  product  by  economic  activities  for  public  sector  at  current 
prices  for  the  years  (2005-2015)  (Million  I.D.). 


Table  5:  Shows  Bayes  Factor  For  H0\  F°0°  4-  co  varsus  Ht:  F0  4-  a) 


Models 

*01 

Evidence 

Model  of  section  1 

1.875279497000141  X  10-25 

very  strongly  favors  H1 

Model  of  section2 

7.903255738497785  x  10“9 

very  strongly  favors  H1 

Total  model 

8.946564506647714  x  10“34 

very  strongly  favors  H1 

Confounded  model 

6.115669635463685  x  10“32 

very  strongly  favors  H1 
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posterior  bO 


posterior  bl 


Figure  1:  Posterior  Density  of  the  Coefficients  for  the  Model  of  Sectionl 
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Figure  2:  Shows  (7000)  Iteration  of  the  Gibbs  Sampler  for  the  Model  of  Sectionl 


Figure  3:  Shows  the  Density  Based  on  (7000)  Iteration  of  o\,o\  and  g\ 
the  Gibbs  Sampler  for  the  Model  of  Sectionl 


Posterior  bl 


Figure  4:  Posterior  Density  of  the  Coefficients  for  the  Model  of  Section2 
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Figure  5:  Shows  (7000)  Iteration  of  the  Gibbs  Sampler  for  the  Model  of  Section2 


Figure  6:  Shows  the  Density  Based  on  (7000)  Iteration  of  o\,a\  and  o\  the  Gibbs  Sampler 

for  the  Model  of  Section2 


Figure  7:  Posterior  Density  of  the  Coefficients  for  total  Model 
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Figure  8:  Show  (5000)  Iteration  of  the  Gibbs  Sampler  for  the  Total  Model 
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Figure  9:  Shows  the  Density  Based  on  (5000)  Iteration  of  a\,a\  And  a\  the  Gibbs  Sampler 


Figure  10:  Posterior  Density  of  the  Coefficients  for  Confounded  Model 
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Figure  11:  Show  (8000)  Iteration  of  the  Gibbs  Sampler  for  the  Confounded  Model 


— I - 1 - 1 - 1 — 

2000  4000  6000  8000 

Iterations 


o 

-2  H 


— I - 1 - 1 - 1 — 

2000  4000  6000  8000 

Iterations 


CM 


i - 1 - 1 - r 

2000  4000  8000  B000 

Iterations 


Figure  12:  Shows  the  Density  Based  on  (8000)  Iteration  of  \  and  g\  the  Gibbs  Sampler 


CONCLUSIONS 

The  conclusions  which  are  obtained  throughout  this  paper  are  given  as  follows: 

•  The  models  for  building  and  construction  and  wholesale,  retail  trade,  hotels  and  others  where  significance  at  a 
0.05  level  of  significance.  This  means  there  is  a  significant  effect  for  gross  fixed  capital  formation  at  current 
prices  for  the  years  (2005-2015)  (Million  I.D.)  for  previous  year  and  gross  domestic  product  by  economic 
activities  for  public  sector  at  current  prices  for  the  years  (2005-2015)  (Million  I.D.)  on  the  gross  fixed  capital 
formation  at  current  prices  for  the  years  (2005-2015)  (Million  I.D.)  for  all  economic  activities 

•  The  total  model  estimator  for  all  the  economic  activities,  is  significance  at  a  0.05  level  of  significance  according 
to  F-  value  which  (20.27009)  with  p-  value  (F)  which  (7.00e-19),  as  well  as  the  value  of  determinant  coefficient 
( R 2  =  0.70),  this  means  the  total  model  is  agree  with  the  economic  theory. 

•  The  total  model  after  confounding  also  significant  at  a  0.05  level  of  significance. 

•  The  values  of  the  parameters  for  the  models  estimates  based  on  the  Bayesian  method,  were  encouraging. 

•  The  Bayes  factors  favors  with  strong  evidence  with  all  models  for  the  data  of  gross  fixed  capital  formation  and 
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